Directed polymers at finite temperatures in 1+1 and 2+1 dimensions 



Xiao-Hong Wang 1,2 , Shlomo Havlin 1 * and Moshe Schwartz 3 
1 Minerva Center and Department of Physics, Bar-Ban University, Ramat Gan 52900, Israel 
2 Department of Thermal Science and Energy Engineering, University of Science and Technology of China, Hefei, Anhui 

230026, P. R. China 

3 School of Physics and Astronomy, Raymond and Beversly Sackler Faculty of Exact Science, Tel-Aviv University, Tel-Aviv 
^\ • University, Tel-Aviv 69978, Israel 

GN ' (February 1, 2008) 

o\ : 

T— I . 

■ We present systematic numerical simulations for directed polymers at finite temperatures in 1+1 

' and 2+1 dimensions. The transverse fluctuations and free energy fluctuations tend to the strong 

/"•n , coupling limit at any temperature in both 1+1 and 2+1 dimensions for long time t. Two different 

definitions for energy fluctuations at finite temperatures, which are the ensemble energy fluctuations 
ON ' and the internal energy fluctuations, are investigated. Apart from zero temperature, the behavior 

of the energy fluctuations and the free energy fluctuations for directed polymers is shown to be 
different. At finite temperatures, the ensemble energy fluctuations in both 1+1 and 2+1 dimensions 
and internal energy fluctuations in 1+1 dimensions scale as i 1 / 2 where the free energy fluctuations 
in 1+1 dimensions and 2+1 dimensions scale as t ' and t respectively. As a consequence of that 
C/j ' the specific heat in both 1+1 and 2+1 dimensions scales as t and the entropy fluctuations in 1+1 

! t ' dimensions scale as t 1 ' 2 at any finite temperature. 

c3 : 

g . Introduction 

' Rough surfaces and interfaces are found in a wide variety of natural and industrial processes. The KPZ equation 
provides a good quantitative understanding of surface growth phenomena 1-5 . It is well known that the KPZ equation 
Q ■ is equivalent to the problem of directed polymers in a random medium 6-9 , which has attracted much attention in 
recent years. The problem of directed polymers is also relevant to many other fields ranging from paper rupture 10 and 
spin glasses 11 to bigh-T c superconductors 12 . The KPZ equation provides the exact dynamical exponent for directed 

' polymers in 1+1 dimesions 6-8 . A self-consistent expansion of the correlation function 13 yields results compatible with 

, simulations for 2+1 dimensions. Above 2+1 dimensions, the behavior of directed polymers has a phase transition when 
t— I ' temperature is raised and the system goes over from the strong coupling behavior to a perturbative behavior 13-18 . 
The strong coupling behavior has been studied by extensive numerical investigations for the energy fluctuation at 
zero temperature 6 ' 8,19-21 , which is ensemble fluctuations of the ground state energy of directed polymers. At finite 
temperatures, fluctuations of the free energy have been studied extensively but it seems that quantities like fluctuations 

• of the internal energy have not received much attentions. 

| In this article, we present systematic numerical simulations for directed polymers in 1+1 and 2+1 dimensions. Our 
numerical simulations show that the transverse fluctuation and free energy fluctuation crossover to the strong coupling 
behavior at any finite temperature in both 1+1 and 2+1 dimensions, in agreement with the predictions of the field 
theory 13,22 ' 23 . We define the energy fluctuations at finite temperatures using two different definitions: the ensemble 
energy fluctuations and the internal energy fluctuations. The difference between the two definitions is related to the 
specific heat. We find that the behavior of the energy fluctuation at finite temperatures is different from that at zero 
temperature. For understanding the energy fluctuations at finite temperatures, we introduce the entropy and study 
its fluctuation. The behavior of the energy fluctuations at finite temperatures can be understood after investigating 
the specific heat and the entropy fluctuations. 



o 



Model 



X, 

Consider a directed polymer on a hyperpyramid lattice structure with the random energy assigned on each bond. 
The partition function G(x,t) for directed polymers starting from (0,0) and ending at (x,t) is defined by G(x,t) = 

^2 er Ec ' l T where Eq is the sum of the energies on the path C and T is the temperature. The iteration relation for 

c 

the partition function G(x. t) is 

G(x, t + l) = G{x- 1, t)e' tl/T + G{x + 1, i)e -e " /T , (1) 
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in which, e; and e r are the energy assigned to the left and right bonds of the point (x,t). The free energy F(t) is 
given by F(t) = —TlnG(f), where G(t) = Y] G(x, t) is the total partition function. The transverse fluctuation Ax 

X 

for the directed polymers is ((x 2 )) 1 / 2 where (x 2 ) — 'J2x 2 G(x,t)/G(t) (A is the ensemble average of the quantity 

X 

A). The transverse fluctuation Ax and free energy fluctuation AF = (F 2 — F )*/ 2 have been commonly studied. 
Another interesting quantity for directed polymers is the energy fluctuation, ft is clear that the energy fluctuation 
at zero temperature is just that of the free energy. For finite temperatures, we propose two definitions for the energy 
fluctuations. One is the ensemble energy fluctuation and the other is internal energy fluctuation. The ensemble energy 

fluctuation is defined by (AE)b = ((E 2 ) — (E) ) 1//2 , in which (E) is the internal energy: 



ljp\ X C rp2^ n G 

= EE e ~ Sc/T = ~9f~ ( ' 

x c 



and 



EEE 2 c e- E o/T 



d r ~dlnG^ r ~dlnG^ 



EE e ~ Bc/T dT <9T dT 



The internal energy fluctuation is defined by {AE)t = ((E) 2 — (E) ) 1 / 2 . It is obvious that both the two energy 
fluctuations will degenerate into the fluctuation of the sum of the energy on the optimal path when temperature tends 
to zero. From the equations (2) and (3), we have 

(AE) 2 B = (AE) 2 T +T 2 C(T,t) (4) 

where C(T,t) = ^[T 2 ^-] = ^ is the specific heat 24 . 

In order to understand the internal energy fluctuation , we use the entropy as S = ((E) — F)/T and the entropy 

fluctuation is AS = (S 2 — S ) 1/<2 . It is easy to see that 

(AS) 2 = ^[(AE) 2 T + (AF) 2 ] - ^B(T,t), (5) 



-2, 



where B(T,t) = (E)F - (E) F = -Ir 3 ^[(AF) 2 /T 2 ]. 

The above quantities can be calculated by the iteration relations given in Appendix. In the following, we present 
systematic numerical simulations based upon the equation (1) and equation (Al) and (A2) from the Appendix with 
the initial conditions G(x, 0) = S x , E(x 1 0) = and E 2 (x, 0) = 0. The random energy assigned on the bond is assumed 
to be uniformly distributed in the interval (-0.5,0.5) and uncorrelated in space and time. The effective program RAN2 
was used for generating the uniform distributed random number 25 . We use length up to t = 2000 (d = 1 + 1) and 
t = 150 (d = 2 + 1). Six thousands configurations for 1+1 dimensions and four thousands configurations for 2+1 
dimensions were collected for the ensemble average. In order to avoid unstability, we normalize the partition function 
at every time step. 

Directed polymers in 1+1 and 2+1 dimensions 

At infinite temperature, directed polymers become random walks: the transverse fluctuation Ax will scales as 
i 1 / 2 and the ensemble energy fluctuation will be (AE) B cx i 1 / 2 because the ensemble energy fluctuation at N time 
steps is only the variance of the sum of N uncorrelated random numbers at the limit. The free energy fluctuation 
for infinite temperature can be predicted by the Edwards- Wilkinson equation 26-28 that describes the linear surface 
growth phenomena. It follows that the free energy fluctuation AF scales as t 1 ' 4 in 1+1 dimensions and tends to the 
logarithmic behavior in 2+1 dimensions. 

When approaching zero temperature, the system will take the strong coupling behavior. Our simulations show 
that the transverse fluctuation scales as t v and both the free energy and energy fluctuations scales as t w , in which 
v = 2/3, u = 1/3 for d = 1 + 1 and v w 0.6, u) w 0.2 for d = 2 + I. These results are in agreement with the theoretical 
treatments and the previous numerical works 6-9 ' 13 ' 19-20 . 

The problem of directed polymers in 1+1 and 2+1 dimensions at finite temperatures is investigated. Our simulations 
show that the transverse fluctuation Ax crossovers from i 1 / 2 to i 2//3 and the free energy fluctuation AF crossovers 
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from t 1 / 4 to i 1 / 3 in 1+1 dimensions at any finite temperature (see Fig. 1) in agreement with earlier arguments 29 . The 
scaling relations Ax = T 2 fi(T~ A t) and AF = Tgi(T~ 4 t), which describe the crossover from the Edwards- Wilkinson 
behavior to KPZ behavior in 1+1 dimensions given by Amar and Family using the scaling analysis 29 , are confirmed 
for T > 0.25 (see Fig. 2). 

Similar to 1+1 dimensions, the transverse fluctuation Ax and free energy fluctuation AF crossover from t 1 / 2 to 
t 6 and from the logarithmic behavior to t 2 respectively in 2+1 dimensions (see Fig. 3). For 2+1 dimensions, we 
can fit the data very well using (AF) 2 = 0.0068Zn(l + 20f) at very high temperature. Assuming the scaling relations 
Ax = Tf 2 (T- 2 t) for T > 0.1 and AF = Tg 2 (T- 2 ln(l + 20*)) for T > 0.2, the data collapse very well (see Fig. 4). We 
note that an argument for the crossover time lnt x ~ T 2 for the free energy fluctuation has been given earlier 16 ' 22-23 . 
Also note that our numerical data suggest that the crossover time for Ax and AF are different. 

The difference between the two different definitions for the energy fluctuation is proportional to the specific heat as 
shown in equation (4). Our simulations shows that that specific heat C(T,t) — > for T — ► as 0(T a ) where a = 1.3 
for d = 1 + 1 ,a = 1.6 for d = 2 + 1 and C(T,t) -> a 2 /T 2 for T -> oo where a = (^) 1/2 is the variance of the 
random numbers. We also find that the specific heat scales as * at any temperature in both 1+1 and 2+1 dimensions 
(see Fig. 5). Fig. 6 shows that the ensemble energy fluctuations (AE) 2 B scale as * as the specific heat apart from 
zero temperature. We find very good fit with the scaling relations (AE) 2 B — e 1 / T hi(e~ 1 ' 5 / T i) in 1+1 dimensions and 
(AF)| = e°- 2/T /i2(e~ - 5/T t) in 2+1 dimensions for T < | (see Fig. 7). Howevever, in both 1+1 and 2+1 dimensions, 
the picture for the internal energy fluctuation is not clear (see Fig. 8). We note that the internal energy fluctuations 
is not monotonic with the temperature. But, at the two limits of zero and infinite temperatures, the internal energy 
and free energy fluctuations are the same in both 1+1 and 2+1 dimensions. For investigating the difference between 
the internal energy fluctuations and the free energy fluctuations, we study the entropy fluctuations. 

Entropy is an important physical quantity for characterizing the disorderd systems. Fig. 9 shows that the entropy 
fluctuation has the behavior (AS) 2 oc * at any temperature in 1+1 dimensions. This result is in agreement with 
earlier arguments of Fisher and Huse 30 . The entropy fluctuation (AS) 2 reaches a maximum about T = T P = 0.2 
and tends to zero as 0(T~ 4 ) when T — > oo and as 0(T) when T — > (see Fig. 10). It follows that the internal 
energy fluctuation and the free energy fluctuation are the same at infinite temperature and the energy fluctuation 
and the free energy fluctuation are the same at zero temperature. The function B(T,t) in the equation (5) can also 
be calculated numerically. From the equation (5), it is easy to see that the function B(T,t) is the same as the free 
energy fluctuation or the internal energy fluctuation at the two limits of zero and infinite temperature. At finite 
temperatures, Fig. 11a shows that the behavior of B(T,t) is similar to that of the free energy fluctuation (AF) 2 and 
also has a crossover from t 1 / 2 to t 2 / 3 in 1+1 dimensions. From the equation (5), we can expect that the internal 
energy fluctuations (AE) T will scale as to t 1 / 2 as the behavior of the entropy fluctuations AS at finite temperatures 
in 1+1 dimensions. 

Discussion 

Using the iteration relations (1), (Al) and (A2), we study the problem of directed polymers in 1+1 and 2+1 
dimensions. The transverse fluctuations and free energy fluctuations have a crossover from the EW behavior to 
KPZ behavior at finite temperatures in both 1+1 and 2+1 dimensions as predicted by the field theory. In order to 
understand the energy fluctuations at finite temperatures, we investigate the specific heat and the entropy fluctuations. 
The specific heat in both 1+1 and 2+1 dimensions and the entropy fluctuations (AS*) 2 in 1+1 dimensions scale as * 
at any finite temperature. As a result, the ensemble energy fluctuations in 1+1 and 2+1 dimensions and the internal 
energy fluctuations in 1+1 dimensions scale as t 1 / 2 at finite temperatures. It means that the energy fluctuations and 
the free energy fluctuations may have a different behavior at finite temperatures in both 1+1 and 2+1 dimensions. 
The behavior of the entropy fluctuations in 2+1 dimensions is more complicated than that in 1+1 dimensions. There 
is an indication of a phase transition in 2+1 dimensions. The phase transition occurs at the temperature where the 
entropy fluctuation reaches the maximum. This result will be published elsewhere 31 . Fig. lib shows that the function 
B(T,t) in 2+1 dimensions also has a crossover from the logarithmic behavior to t 2 at any finite temperature. As a 
result, the internal energy fluctuations in 2+1 dimensions also have a phase transition even if Fig. 8b can not show 
it clearly. 
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Appendix: Iteration relations 

As in deriving the iteration relation for the partition function G(x,t) given by Eq. (1), we concentrate on the 
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1+1 dimensional case. The generlization to higher dimensions is obvious. In order to obtain the iteration re- 
lations for calculating the energy and entropy fluctuations, we define E(x,t) = E C ( x ^e~ Ec ^^/ T / G{t) and 

C{x,t) 

E 2 (x,t) = E l( x t) e~ Ec ^^ T /G(t). It is clear that (E) = E-^OM) and (E 2 ) = J2E 2 (x,t). The derivation for 

C(x,t) ' x x 

the iteration relations for E(x, t) and E 2 (x, t) are given as follows: 
E(x,t + 1)= £ %,i+i)e" Bo( " H,)/T /G(*+l) 

C(x.t+1) 

= [ E (Ec(x-i.t) + ei)e-^-^)/ T + £ (^c(x+i,t) + e r )e-^+^>+^/ T ]/G(t + 1) 

C(x-l,t) C(x+l,t) 

= [e-^/ T E(x - 1, t)G(t) + e- e -/ T ^(x + 1, t)G(t) + e ;e - £ </ T G(x - 1, t) + e r e~^' T G{x + 1, t)]/G(t + 1) (Al) 
and 

1) = E 4( llf+ i)«" W MHl) 

C(x.t+1) 

= [ E (^c(x-i,t)+Q) 2 e- (£c '(- 1 - t ) +£,)/r + E {E c(x+ht) +e r fe-^ + ^r)/T ]/G{t+1) 

C(x-l,t) C(x+l,t) 

= [e-^' T E 2 {x - 1, t)G(i) + e-^/ T i? 2 (a; + 1, t)G{t) + 2e /e -^/ T i?(a; - 1, t)G{t) + 2e r e-^/ T ^(a; + 1, t)G(t) 

+efe-^/ T G{x - 1, t) + 4e-^ T G(x + 1, t)]/G(t + 1). (A2) 
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FIG. 1. (a) Plot of the transverse fluctuations (A:r) 2 as a function of time t at different temperatures T = 10, 1, ^, |, ^ 
(from bottom to top) in d = 1 + 1 dimensions, (b) Plot of the free energy fluctuations (AF) 2 as a function of time t at different 
temperatures T = 10, 1, |, |, |, (from bottom to top) in d = 1 + 1 dimensions. 




FIG. 2. Scaling plot for the transverse fluctuations (Aa:) 2 shown in Fig. la and the free energy fluctuations (AF) 2 shown in 
Fig. lb for d = 1 + 1 dimensions where circle, square, diamond, triangle up, triangle left, triangle down correspond to different 
temperatures T = 10, 5, 2, 1.0, 0.5, 0.25 respectively. 
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FIG. 3. (a) Plot of the transverse fluctuations (A:r) 2 as a function of time t at different temperatures T = 1, |, i, -jj 
(from bottom to top) in d = 2+ 1 dimensions, (b) Plot of the free energy fluctuations (AF) 2 as a function of time t at different 
temperatures T = (from bottom to top) in d — 2 + 1 dimensions. 
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FIG. 4. (a) Scaling plot for the transverse fluctuations (Aa;) 2 in d = 2 + 1 dimension where circle, square, diamond, 
triangle up, triangle left, correspond to different temperatures T — 1, |, i, i, respectively, (b) Scaling plot for the free 
energy fluctuations (AF) 2 in d = 2 + 1 dimensions where circle, square, diamond, triangle up, triangle down, star, triangle 
right correspond to different temperatures T = 5, 1, |, |, i, i respectively. 
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FIG. 5. Plot of the specific heat as a function of temperature, (a) for the different time t — 100, 500, 1000, 2000 in d = 1 + 1 
dimensions and (b) for the different time t — 50, 100, 120, 150 in d — 2 + 1 dimensions. 




FIG. 6. Plot of the ensemble energy fluctuations (AE) B as a function of time t, (a) for different temperatures 
T = 10, |, |, |, |, |, ^ (from top to bottom) in d — 1 + 1 dimensions and (b) for different temperatures T — l,|,i,i,-i,^,i 
(from top to bottom) in d = 2 + 1 dimensions. 
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FIG. 7. (a) Scaling plot for the ensemble energy fluctuations (AE)% in d = 1 + 1 dimensions where circle, square, diamond, 
triangle up, triangle left, triangle down correspond to different temperatures T = -L, |, i, i, | respectively, (b) Scaling plot 
in d = 2+ 1 dimensions where circle, square, diamond, triangle up, triangle left correspond to T = fg, jg, -j^, § respectively. 
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FIG. 9. Plot of entropy fluctuations (AS) 2 as a function of time t in d = 1 + 1, (a) for different temperatures 
T = |, |, yjj, jg, -L, -L, ^j, i (from top to bottom) and (b) for different temperatures T = 10,5,2, 1, |, |, |, | (from bot- 
tom to top). 




FIG. 10. Plot of entropy fluctuations per unit length of the polymer (AS) 2 /t as a function of temperature T for the different 
time t = 50, 500, 1000, 2000 in d = 1 + 1 dimensions. 
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FIG. 11. Plot of the function B(T, t) given in the equation (5) of the text as a function of time t, (a) for different temperatures 
T = 10, 1, |, |, |, |, ^ (from bottom to top) in d = 1 + 1 dimensions and (b) for different temperatures T = 1, |, |, i, i, |, ^ 
(from bottom to top) in d = 2 + 1 dimensions. 
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